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ABSTRACT 

Aims. To demonstrate the capabilities of regularized inversion to recover Differential Emission Measures (DEMs) from multi- 
wavelength observations provided by telescopes such as Hinode and SDO. 

Methods. We develop and apply an enhanced regularization algorithm, used in RHESSI X-ray spectral analysis, to constrain the 
ill-posed inverse problem that is determining the DEM from solar observations. We demonstrate this computationally fast tech- 
nique applied to a range of DEM models simulating broadband imaging data from SDO/AIA and high resolution line spectra from 
Hinode/EIS, as well as actual active region observations with Hinode/EIS and XRT. As this regularization method naturally provides 
both vertical and horizontal (temperature resolution) error bars we are able to test the role of uncertainties in the data and response 
functions. 

Results. The regularization method is able to successfully recover the DEM from simulated data of a variety of model DEMs (single 
Gaussian, multiple Gaussians and CHIANTI DEM models). It is able to do this, at best, to over four orders of magnitude in DEM 
space but typically over two orders of magnitude from peak emission. The combination of horizontal and vertical error bars and the 
regularized solution matrix allows us to easily determine the accuracy and robustness of the regularized DEM. We find that the typical 
range for the horizontal errors is AlogT w 0.1 - 0.5 and this is dependent on the observed signal to noise, uncertainty in the response 
functions as well as the source model and temperature. With Hinode/EIS an uncertainty of 20% greatly broadens the regularized 
DEMs for both Gaussian and CHIANTI models although information about the underlying DEMs is still recoverable. When applied 
to real active region observations with Hinode/EIS and XRT the regularization method is able to recover a DEM similar to that found 
via a MCMC method but in considerably less computational time. 

Conclusions. Regularized inversion quickly determines the DEM from solar observations and provides reliable error estimates (both 
horizontal and vertical) which allows the temperature spread of coronal plasma to be robustly quantified. 

Key words. SumCorona - Sun:Flares - Sun: X-rays, gamma rays - Sumactivity - Sun:UV radiation 



1. Introduction 

Observations of the solar atmosphere with temperature sensi- 
tive spectral lines provide crucial information about the tem- 
perature distribution of the emitting plasma. These are vital 
for trying to resolve the question of which mechanisms heat 
different solar phenomena. Such as whether coronal loops are 
heated by a nanoflare mod el of magneti cally reconnecting multi- 
braided loo p strands (e.g. | Parkerf [l988) or chromospheric evap- 
oration (e.g. lHiravamall974l) . Or is the hot emission observed in 
large through to micro- flares dFletcher et al.l20lHlHannah et alj 
|201 lb predominantly due to energetic particles or other mecha- 
nisms such as waves. To reliably answer these questions, one 
needs to know not only the uncertainties on the emission for a 
given temperature, but also the uncertainties on the temperature 
itself, i.e. the temperature resolution. 

T he observations made with Hinode's X-ray Telescope 
XRT riGolubet al]l2007l) and EUV Imaging Spectrometer EIS 
dCulhane et al.l I2007I) an d SDO's Atmospheric Imaging Array 
AIA dLemen et al|201 lb and EUV Variability Experiment EVE 
(I Woods et al.H2010h provide a wealth of information about the 
solar emission over a broad range of temperatures. Assuming 
this UV/EUV/X-ray emission is both optically thin and in ther- 
mal equilibrium, via collisions, then the temperature distribu- 



tion of plasma emitting along the line of sight h can be de- 
scribed by the differential emission measure DEM, typically 
given by %(T) = n 2 dh/dT [cm~ 5 K _1 ] where n(h(T)) is the elec- 
tron de nsity at h and with temperature T (see Chapter 4 Mariska 
(1992) for detailed discussion of the different DEM forms). 
This however, can not be immediately inferred from such multi- 
wavelength observations as the DEM is convolved by the emis- 
sion processes and the instrumental response, i.e. 
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where gj is our observable for the 2 th filter, which has a tem- 
perature dependent response function Kff) and 6gi is the er- 
ror. For spectroscopic observations this is respectively the line 
intensity and contribution function (examples shown in Figure 
[TJ. The uncertainties associated with these observations (count- 
ing statistics, background and instrumental errors) compounds 
the difficulty in the determ ination of the DEM and results in an 
ill-posed inverse problem dTikhonoyl 1963t iBertero et al.l tl985: 
ICraig & Brownlll986t ISchmittet al.lll996l) . Any direct attempts 
to solve Eq. Q] normally leads to the amplification of the uncer- 
tainties, and hence spurious solutions. 

To reconstruct the DEM additional information (i.e. con- 
straints) has to be added and numerous approaches have been 
developed to solve this problem (for overviews see, for instance, 
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Monsigno ri Fossi & Landinl (Il99ll) or Chapter 5 Phillins 'et alJ 
(I2008h 1. The simplest of these is to assume that all the emission 
is at a single temperature (isothermal) with £(T) oc S(T - To) 
where 6(x) is the Dirac delta function. The ratio of emission 
between two filters is then equal to the rat io of the response 
functions at the i sothermal temperature (e.g. IWeber et al.| [2005; 
Real e et al.l2009l) . Dividing the observable by the response func- 
tion gi/Kt and plotting this as a function of temperature, the 
intersect point of the different curves (EM loci curves) will 
give the isothermal temperature and emission measure (e.g. 
ISchmelz et alj20l1b . Although this is a simple and computation- 
ally fast method it does require the isothermal assumption and if 
the DEM is multi-thermal this method will produce erroneous 
results. 

Another approach is to forward fit a chosen model, 
minimising the differences in observable space. This has 
been implemented for a d is cretise d spline mode l DEM 
(iMonsignori-Fossi & La ndin 3 119921: iBrosius et al.l [l996; 
iParenti et al.l |2000|) and more recently using the IDL mpfit 
routine from SDO/AIA and Hinode/XRT in M. Weber's 
xrt_dem_iterative2 .prcQ. dWeber et al.l l2004t iGolub et all 
This iterative forward fitting approach has also been 
developed with mu ltiple Gaussian model D EMs using the IDL 
POWELL routine (lAschwanden & Boernerll201 lh . To estimate 
the error in the DEM with these methods a Monte Carlo 
approach is adopted, producing multiple realisations within a 
given noise range. These approaches will find parameters for 
the model DEM but requires an assumed model and can be 
computationally slow when error estimates are required. The 
Metropolis MCMC approach has also been used on recovering 
the D EM for a large set of EUV spectral lines dKashvap & D rake 
1998), part of PINTofALE spectral analysis package. This will 
give a robust measure of the parameter probability space but 
again can be computationally intensive, especially if a large 
number of lines and model parameters are considered. Bayesian 
formalism has also been recently used in a Bayesian Iterative 
Method (BIM), successfully reconstructing DE Ms from both 
simulated and observed data ( Gorvaev et al. I20T0I) . The iterative 
mpfit method was compared to a maximum likelihood and a 
genetic algorit hm technique, findin g similar results between the 
approaches (Siarko wski et alJ[2008l) Another genetic algorithm 
approach, which involves a preconditioning step where the 
optimum subset o f spectral lines are s elected, was found to 
be more effective (Mcintosh et al. 2000). Currently a Singular 
Value Decomposition (SVD) inversion approach is also under 
development by Weber for SDO/AIA data. 

Regularized inversion methods introduce an additional 
"smoothness" to constrain the amplification of the uncertainties, 
allowing a s table inversion to reco ver the DEM solution (e.g. 
ICraidll977t ICraig & Brownl ll986). This was demonstrated to 
have promise for solar observation s bvlCraigl d!977|) and subse- 
quently tested on simulated data by Monsignori Fossi & Landini 
d!991h . Several forms of regularized inversion - truncated (or 
"zeroth-order") SVD, second-order regularization and maxi- 
mum entropy regularization - hav e also been tested u sing simu- 
lated EUV spectral line emission dJudge et alJl!997h . Although 
they determined that these approaches were superior to the sim- 
ple ratio method they found several problems with the regular- 
ized inversion: the smoothness criterion used may not be phys- 
ically "appropriate"; the solutions are highly sensitive to uncer- 




5.0 5.5 6.0 6.5 7.0 7.5 

Log, (T) [K] 




5.5 6.0 6.5 

Log, (T) [K] 



Fig. 1. Temperature response for SDO/AIA ( Boerner et al l20lTb 
and the contribution functions for several EUV lines fr om 
CHIANTI observable by Hinode/EIS dLandi & Young|2009l) . 



tainties in the kernel (K, response or contribution functions); the 
return of negative solutions. 

In this paper we present a regularization method, which re- 
solves some of these problems and robustly recovers the under- 
lying DEM with errorfl The method not only determines the 
DEM quickly, and its associated errors, but also naturally pro- 
vides an estimate to the temperature resolution of the method. 
The ability to quickly compute (via Generalised Singular Value 
Decomposition) the DEM and its associated errors is due to 
the fact this method is linear, unlike those using maximum 
entropy for example. This method has already been imple- 
mented and ap plied to solar data for the inversion of RHESSI 
dLin et alJ 2002f) X-rav spectra to their source electron distribu- 
tion ( Kontar et al.l20 04). Several other inversion techniques have 
been developed to infer X-r ay photon spectra and/or electron 
spect ra from RHESSI data dPiana etal J 120031; iMassone etakl 
2004; iBrown et alf 2006) however it is the regularized inversion 
that has become the de facto approach. It has subsequently been 
used to infer the DEM, as wel l as the non-the rmal emission, from 
RHESSI hard X-ray spectra dPrato et alJl2006l) . 

In §2 we detail the regularization method and how the er- 
rors and temperature resolution are determined. In §3 and §4 
we demonstrate the capabilities of the regularization method, 
in comparison to other methods, on simulated SDO/AIA and 
Hinode/EIS data for a variety of model DEMs (single and multi- 
ple Gaussians and the CHIANTI DEM models). In §5 we use the 
regularization method to recover the DEMs from observations of 



1 Available in SolarSoftWare with the Hinode/XRT software 
$SSW/hinode/xrt/idl/util/xrt_dem_iterative2.pro 



2 The code written in IDL requiring 
http://www.astro.gla.ac.uk/~iain/demreg/ 



SSW is available online: 
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i log T[K] log T[K] 

Fig. 2. Singular values 0, (left), singular vectors w, (middle) and the regularized solution Eq.[6]as a function of i (right). The lines 
in the singular value plot (left) indicate the regularization parameter A that produces the solution with desired x 2 (dashed line) and 
with the additional positivity constraint (dash dot line). These values are from the regularized solution shown in Figure|3] 



an ac tive region made with Hinode/EIS and XRT dWarren et alJ 
2010). 



2. Regularized inversion of multi-wavelength data 



To find the line of sight DEM £ (Tj) j 
solve the system of linear equations: 



l,..,M[cnr 5 K-']is to 



(2) 



where gi is the observable (either the imaged DN or the in- 
tegrated line intensity) for the specific filter or wavelength i 
(i — 1, N) and K,-. ; is the corresponding temperature response 
or spectral line contribution function (examples shown in Figure 
[TJ- Eq.[2]is a generally ill-posed inverse problem and hence the 
least square problem 



Kf(r)-g 



8g 



(3) 



does not have a unique solution. In the case M = N, a formal 
solution of Eq.|2]can always be written as f(7j) = Krjg,-, where 

K -1 is the inverse of matrix K. However, due to the large con- 
dition number, which is the ratio of the largest to the smallest 
singular values of K, such a solution is meaningles s in all prac- 
tical cases due to substantial noise amplification (ICraipl fl977t 
iBertero etal.1119851) . In addition, with M > N the system is also 
under-determined. With inherent statistical and/or instrumental 
noise/uncertainties in the data g, the information about the true 
solution %(T) is lost and cannot be recovered without adding ex- 
tra information about %(T). Therefore constraints must be ap- 
plied to obtain a unique meaningful solution. All methods solv- 
ing this system to find the DEM %(T) explicitly or implicitly add 
information not present in the data to obtain the approximate so- 
lution. The simplest, but most popular way to constraint the data 
is to fit a model function %(T, <?,) with a number of free parame- 
ters a, that minimise the Eq.[3] Forward fitting is highly unsatis- 
factory if the functional shape of %(T) is a priori unknown. 

As any attempt to reconstruct the DEM directly leads to 
substantial noise amplification in £(T), the broad approach to 
achie ve a solution is to add linear cons t raints to the DEM 
(e.g. lTikhonov|ll963t IBertero et all I1985L 11988b ICraigl Il977t 
ICraig & Brownl 19861) Often, so-called zero order regularization 
is used, which selects the smallest norm solution out of infinitely 



many possible solutions. This approach pro ved to be robust for 
vario us problems and is not over restricting dBertero et all 1985. 
[19881) . Hence, we solve the least square problem 

||Kf(r) - gf = min subject to ||L(f(T) - fr(T))\\ 2 < const, (4) 

with K = (£g)~'K and g = (5g)~ l g. This can be solved using 
Lagrangian multipliers, i.e. 



||K f(D - -gf + A\\L (f(D - MT))\\ 2 = min, 



(5) 



where L is the constraint matrix, A is the regularization param- 
eter (related to the^ 2 of the solution), and %o(T) is the "guess" 
solution, which will be explained in detail below. The L2 norm 
is defined as a sum ||x|| 2 = x r x = Yli=\^ over a U filters or 
intensities. Importantly, the solution of Eq.[5]is unique and well- 
behaved. The formal solution of Eq. |5]^(r) can be simply ex- 
pressed in matrix form as a function of regularization parameter 
A but to avoid time consuming matrix ma nipulations G eneralized 
Singular Value Decomposition GSVD dHansenlll992h is used. 
The GSVD of matrices K e R MxN , L 6 R NxN produces a 
set of singular values y,, /?, and singular vectors u,, v,, w,, with 
i = l,...,N which satisfy y] + /3 2 = 1, U T KW = diag(y) 
and V T LW = diag(j6). These then provide the solution (Hansen 
1 19921) to the minimization problem given in Eq.[5]as 



f(g-u,) W/ , ^b(D 



tf + A 



Ji 



n 



(6) 



with <pi = yilPi- This solution weights the contributions from 
various singular vectors differently, filtering out the singular vec- 
tors with i, for which (f> 2 < A. Henc e removing un-phy sical os- 
cillatory component of the solution dBertero et al.ll 1 988b . Figure 
|2]shows typical behaviour of singular values and vectors as well 
as the construction of the solution. 



2.1. Regularization Parameter 

To find the solution %a(T) is to determine the regularization pa- 
rameter A, which is done using Morozov's discrepancy principle 
dMorozovll967l) . i.e. 



jj ||K £i(r) 



■8 



(7) 
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solving for A after substituting £a(T) from Eq. [6] Here a is 
the regularization "tweak" value which effectively controls the 
required^- 2 of the solution in observable space. The a value has 

a clear meaning when ||K £a(T) - gll are normally distributed, 
with a mean of N and variance of y/N. Therefore values of 
||K £a(T) - gm in the range N - s/N < N + yfN are accept- 
able values and a helps choose the exact value within this range. 
This also helps to put more or less weight on the data, with a < 1 
requiring a "better" agreement with the data. 

2.1.1. Positively defined DEM 

Using the method discussed in the previous section, it is pos- 
sible to select only positive solutions from the family of solu- 
tions %a(T) by choosing an appropriate A. The method based on 
Morozov's discrepancy principle chooses the parameter A that 
gives the DEM solution ijz(T) with the desired a in Eq.|7] The 
intrinsic DEM from the Sun should be positive but the regu- 
larization method provides no guarantee of a positive solution. 
Although this appears to be a problem with this implementa- 
tion of the regularization method, one should remember that the 
DEM derived from observations need not be positive given the 
often poorly known uncertainties, response functions and the 
possibility of background subtraction. 

A positive DEM solution can be achieved with an additional 
criterion to the choice of regularization parameter. That is we 
take the regularization parameter A that provides the smallest 

||K^(r)-g|| 2 - aN and fr(T) > 0. This approach has the 
advantage of maintaining the linear calculation of the solution 
unlike those that try to implement the positivity constraint in 
Eq. Indir ectly, producing non-linear or iter ative solutions (e.g 
iPiana & Bertero1ll997blde Vilhers et al.lll999h . As a is the x 2 of 
the regularized solution, the ability to recover a positive solu- 
tion strongly depends on the error estimates on the input data 
and knowledge of K. If the error used is too small the x 1 of 
a positive solution can be erroneously high and in general, the 
positivity constr aint produces a larger y 2 , behaviour previously 
demonstrated by Bertero & Dovil ([T98 1 ). 

2.2. Initial Guess Solution 

The standard mode of operation of our regularization algorithm 
requires no initial guess solution, i.e. £o(T) = 0. However, %o(T) 
is used in the calculation of the constraint matrix (either for the 
higher order constraints or in the constraint weighting, see 32.31 ). 
To avoid this problem we run our regularization algorithm (solv- 
ing Eq. [7] and then Eq. [6} twice. On the first run the guess so- 
lution is %o(T) = 0, the constraint matrix is the identity matrix 
L = I and we find a weakly regularized solution with a = 10. 
The regularized solution found can then be used as the initial 
guess for the second run (£o(T) - £r(T)) and in calculating the 
chosen constraint matrix (as discussed in 92.3l l. For this second 
run a stronger regularization is used a — 1 to find the final solu- 
tion, then the associated errors in the DEM and temperatures are 
calculated, see 32.41 

An alternative approach can be taken when working with 
high resolution spectroscopic line data as we can use the mini- 
mum of the EM loci curves as the initial guess solution. The EM 
loci curve for each spectral line is the ratio of the line intensity 
to contribution function w gi/K/. As this estimates the EM based 
on the isothermal temperature, it provides the upper limit to the 
DEM as a function of temperature. Multiple spectral lines across 



a wide range of temperatures provides a strong constraint to the 
DEM space and a useful initial guess solution. When selected it 
is automatically calculated within our code and as we are start- 
ing with a non-zero guess solution, this approach only requires 
one run of the algorithm. These two approaches to the guess so- 
lution will be demonstrated for Hinode/EIS data, simulated from 
a variety of model DEMs, in §3.2 and §4.2. For broadband data 
(such as those from SDO/AIA), only the ^o(T) = approach 
is used as the minimum of the EM loci curves provides a poor 
guess solution, especially when the DEM is not isothermal. 

2.3. Constraint matrix 

There is a number of different choices for L and here we con- 
sider only linear constraints. Physically, the quantity J n 2 e dh = 
J £( T)dT is the total number of electrons along the line of sight, 
so_[]£(r)|| 2 < constant, similar to the X-ray case jPiana et al.l 
120031) . This corresponds to the constraint matrix of Lo oc I, a 
zeroth-order constraint. So, applying zero order regularization, 
we find DEM with the smallest amount of plasma required to 
explain the observational data. The source averaged ^(T) results 
from a combination of heating, cooling and and the physics of 
heat transport in the radiating source. Therefore, when %(T) is 
the solution of some differential equation, one can expect that 
g(T) should be differentiable or equivalently the constraint ma- 
trix L is the first order or second order der ivative (Li ~ D 1 o r 
L2 ~ D 2 ), again similar to the X-ray case dKontar et al.l l2004). 
Hence higher order regularization constraints select solutions 
with the smallest variations in temperature and are more restric- 
tive then zeroth-order solutions. The more restrictive methods 
put more weight on the a priori constraints rather than the anal- 
ysed data sets, which could be advantageous for poorly mea- 
sured data. 

2.4. Error and Temperature Resolution of DEM 

Suppose the true DEM solution ^t me {T) is given. Then we can 
write 

g = K£ rae (r) + Sg. (8) 

Any regularized solution to a linear problem can be viewed as 
the replacement of the generalized inverse K + with the regu- 
larized inverse R\ (with liiru^o Ri = K + ), so that our regular- 
ized solution is %r(T) = R^g. Indeed for £o(T) = 0, Eq. [8]can 
straightforwardly re-written in such form with R,i expressed via 
know n GSVD vectors a nd values and exists for all linear meth- 
ods dBertero et al.l[l9 88S). On the other hand, the true solution to 
Eq. [8] can be formally written K + g. Therefore, to estimate the 
error we find the difference between the true solution and our 
solution 

6&T) = &(!■) - U,(T) = (R,K - m me (T) + R,i<5g (9) 

Eq. [9] shows the important result that the error comes from two 
parts: the last term gives us that the noise propagation (vertical 
errors) and the first term gives the temperature resolution (hor- 
izontal errors). This equation presents the method independent 
definition of both the horizontal and vertical resolutions . While 
the noise propagation is normally accounted by the DEM meth- 
ods, the resulting temperature resolution is often not considered. 

The error on the DEM 6%{T) (the vertical error R^^g) is cal- 
culated using the standard Monte Carlo approach ([Press et al.l 
11992^ IPiana et al.1 12003b iKontar etal.1 l2004t iPrato et al.l 12006 j 
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Fig. 3. Results of DEM reconstructions from simulated SDO/AIA data of a single Gaussian model. Shown is the zeroth-order 
regularized solution (left panel, red error bars) and with an additional positivity constraint (middle panel, blue error bars). The 
resulting simulated SDO/AIA data for each filter and the residuals of the regularized solutions are shown in the upper and lower 
right panels. Also plotted is the solution from xrt_dem_iterative2 . pro (grey histogram, left panel) and the EM loci curves 
(purple lines, left panel). 



with multiple random realisations of g within the noise range 
<Sg. Then the one sigma spread of the regularized solutions from 
these realisations provides the measure of the uncertainty on the 
DEM. This is possible due to the linear nature of Eq. [9] In gen- 
eral, the exact statistics of the errors is needed but as this distri- 
bution is unknown, we assume a Gaussian distribution of stan- 
dard deviation 6g for each filter. Then the the vertical uncertainty 
on the DEM are calculated as the standard deviation using 300 
MC Gaussian noise realisations. However, as the probability dis- 
tribution of 6g does not need to be Gaussian, the resulting DEM 
uncertainties could also be non-Gaussian. 

The temperature resolution (horizontal error) of any linear 
inversion method - or temperature bias of the solution - is how 
much the product R^K differs from the identity matrix I. When 
the regularized inverse is similar to inverse of the kernel matrix 
Ri ^ K + , the temperature resolution does not degrade giving 
R^K ^ I. However, for the ill-conditioned problem of DEM de- 
termination, R^K is not identity matrix, but has a finite spread. 
The temperature resolution is then simply the FWHM of R^K 
for a given temperature bin. This represents the temperature bias 
measure or the smallest temperature difference which can be 
meaningfully distinguished in the solution. Conveniently, R^K 
is easy to calculate from the singular values and vectors ob- 
tained in GSVD decomposition of K and L used to find the reg- 
ularized solution. Namely, R^K = WYW 1 where the column 
vector w, forms matrix W and Y is a diagonal matrix with the 
elements constructe d with singular values Y i; = y?/ (y? + A0f) 
dKontar et alj|2004l) . 

When R^K = I the response is diagonal or impulse-like for 
a given temperature, and therefore the true DEM £; tme (T) is not 
distorted. In any practical situation, R,)K is not diagonal but has 
off-diagonal elements (an example is shown in Figure|4]of a sim- 
ulated DEM discussed in detail in §13. U . The off-diagonal terms 
are a spread about a peak value at the diagonal (the red exam- 
ple in Figure @}, so for each temperature Tj the row of (R^K), 
has a finite width. The FWHM of this spread is then taken as 
the temperature resolution for that particular temperature. When 
the off-diagonal terms dominate (the blue and green examples 
in Figure |4]l, the FWHM can still be used to indicate the poor 
temperature resolution. It should be noted that R^K depends on 
the errors 5g, so larger errors results in a poorer temperature res- 
olution, e.g. wider row of the matrix R,iK. The dependency on 
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Fig. 4. (Left) R^K matrix with each row providing the tem- 
perature resolution information for each temperature Tj. The 
highlighted rows are plotted in the right panel also showing the 
FWHM estimate used to calculate the temperature resolution. 
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Fig. 5. Effect of changing the temperature bin size for the single 
Gaussian model DEM, shown in Figure [3] 



6g comes via the regularization parameter A given by Eq.[7] This 
definition of temperature resolution is more conse rvative as it 
does not assume a form of DEM as in IWebeii(l2009l) . 



2.5. Temperature response and instrument uncertainties 

One of complicating factors is that the kernel (the matrix K) of 
the integral equation to be inverted is only known to a limited de- 
gree. The errors come from the calibration uncertainty of the in- 
strument itself and the uncertaint y in the dominant spectral con- 
tribution of each bandpass (e.g. lAschwanden & Boernerfl20l"ll : 
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lO'Dwver et alJEoTll) . Therefore, the temperature response has 
an uncertainty 6K and the linear problem, Eq.[HJ becomes 



g = (K + 6K)^ me (T) + 6g 



(10) 



This translates into an additional uncertainty for £(T) compared 
toEq.|9] 



6£(T) = (R,,K - I)ft rue (r) + R,K 8g + R A 6K ^ e (T) 



(11) 



When the uncertainty (5K is dominated by a systematic error in 
the intensity of the lines measured then it is not temperature de- 
pendent but a constant scaling factor per filer or line (i.e. rows of 
K). Since the shape of the response or contribution function does 
not change as a function of temperature then this is identical to 
the introduction of an additional error to 5g with Eq.fTTIreducing 
to Eq.|9] This will be investigated further in 33.2l and 34. 21 with 
reference to Hinode/EIS. 
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3. Simulated Data: Gaussian Model 

We test the regularization method on simulated SDO/AIA ( {33. Il l 
and Hinode/EIS data ( 33. 21 ) of Gaussian model DEMs. These 
have the form 



No 



y/2ncr T 



exp 



-(log Tj- log To) 2 



2<X 2 



(12) 



where log To is the centroid temperature, cr T is the standard de- 
viation and A^o = f £(T)dT. We consider model DEMs of one to 
three Gaussian components. Note that as we calculate everything 
using log T instead of T, a conversion factor of T In (10)iilog T 
is required in Eqs.Q~]and[2] 

3.1. SDO/AIA Simulated Data 

To simulate SDO/AIA data we take the model DEM from Eq.fHl 
and calculate the expected observable signal in ea ch of AIA's six 
coron al filters g, using the response functions Ki dBoerner et al.l 
12011b . shown in the top panel of Figure Q] From this we cal- 
culate the associated error 6gi using the readout noise and pho- 
ton counting statistics (correcting from DN s -1 px~' to photons 
via the electron and photon gains). Then Gaussian noise within 
these 6gi is added to gy. Our simulated observables g, 5g and the 
response functions K are the inputs to the regularization algo- 
rithm. As described in $2] we run the regularization twice: on 
the first run we use a = 10, L = 1,£q(T) = which provides a 
guess solution that can be used to calculate the desired constraint 
matrix Lo, Li, L2 and final regularized solution with a = 1. 

We first consider a single Gaussian model of A^o = 3.76xl0 22 
cm -5 , cr T = 0.15 and log To = 6.5, shown in Figure[3] Here the 
regularized DEM (red error bars) was found using Lo (zeroth- 
order constraint) and well matches the original model DEM 
(black dashed line). In observable space (right panel in Figure 
the residuals between the model and regularized solution are 
small with^- 2 * 1.0 close to the desired value set by a — 1. Here 
X 2 is taken as the sum of the square of the residuals divided by N 
to match the version in Eq.|7] The R^K matrix for this regulariza- 
tion (Figure[4]i is diagonal over approximately log T — 6.1 -6.9. 
The horizontal error bars (temperature resolution) at each tem- 
perature bin is taken as the standard deviation of the FWHM of 
the rows of R^K. Outside this range the matrix is clearly not 
diagonal, producing large horizontal errors, so the regularized 
DEM is not reliable at these temperatures. 
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Fig. 7. Testing the effect of changing the constraint matrix or- 
der (Lo, Li, L2, zeroth, first and second order respectively). The 
top four panels show the different recovered DEM and result- 
ing residuals for the single Gaussian model shown in Figure [3] 
The bottom four panels shows the same for the broad double 
Gaussian model shown in Figure [6] 



At the lowest temperatures (log T < 6.0) some negative 
DEM values are found, not shown on the log-scale, but this is 
in a temperature range where the DEM contribution is minus- 
cule (over 4 orders of magnitude smaller than the peak) and the 
temperature response is very weak. At the highest temperature 
(logT > 7.0) we have positive regularized DEM but with large 
horizontal and vertical errors again due to being in a range where 
the DEM component is very small. In the middle panel of Figure 
|3]we show the regularized solution (blue error bars) in which the 
regularization parameter has been chosen to minimise x 2 to the 
desired value whilst forcing a positive DEM. The result is a very 
slightly different DEM, particularly below logT w 6.1, and is 
well within the error bounds shown for both regularized solu- 
tions. As expected the^- 2 of the solution is higher (5.5 instead of 
the desired 1) but within the DEM error bounds and so both are 
valid solutions. 

For comparison we have calculated the DEM solution using 
the iterative forward fitting routine xrt_dem_iterative2 .pro 
(black histogram in left panel Figure [3} with 75 Monte Carlo 
(MC) realisations of the solution within the observable error 
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Fig. 6. Example DEM from zeroth-order regularization of simulated SDO/AIA from a two Gaussian model. (Left) DEM space, 
(middle) observable space and residuals (right) R.iK used to calculate the temperature resolution. Also plotted is the solution from 
xrt_dem_iterative2 . pro (grey histogram, left panel) and the EM loci curves (purple lines, left panel). Note that the solution here 
is positive without the additional constraint. 



bound gi + 6gi (grey histograms). This also agrees well with the 
model DEM and also shows a larger spread in DEM solution in 
the temperature ranges where the regularized solution has large 
errors. Again this is due to the minor contribution to the DEM 
(about four orders of magnitude smaller than the peak value) in 
these temperature ranges. Also for comparison we have plotted 
the EM loci curves (gi/(KjTdlnT)) which do not intersect at a 
single point as this is not an isothermal model DEM. The DEM 
solutions are below the EM loci curves, as expected, since these 
estimate the upper limit of the emission. 

The error bars shown for the regularized solution are not in- 
dependent of each other and map out an error boundary region. 
This can be seen in Figure [5] where we show the regularized 
DEM using a variety of temperature binsizes. To achieve this 
we need to interpolate the temperature response functions (top 
panel FigureQ]i from the original binning to our chosen tempera- 
ture binsize. Increasing the number of temperature bins does not 
change the shape of the regularized DEM but produces a clearer 
definition of the error bound, with the overlapping error bars in- 
dicating that the nearby DEM bins are clearly not independent. 
This does however slow down the computation of the regular- 
ization process, although at worst it still only took a few seconds 
to compute. Note again that the vertical errors are taken as the 
variance of 300 regularized solutions found from random reali- 
sations within gi + 5gi and so we are explicitly assuming a simple 
Gaussian spread of vertical error in the regularized solutions. In 
reality the distribution of these errors will be more complicated. 

A second example is shown in Figure [6] where an additional 
broad Gaussian component (A^o = 8.77 x 10 22 crrT 5 , crj = 0.5 
and log To = 6.5) has been added to the model DEM shown in 
Figure[3] Immediately it is clear that the regularization produces 
a RjK that is diagonal over a wider temperature range, indi- 
cating that the narrow temperature range found for the single 
Gaussian model (Figure [3]) was mostly due to the DEM model 
dominating over a small temperature range. For the broader 
DEM model there are still deviations from a diagonal R^K at 
the ends of the temperature range chosen but this is due to the 
limited response of AIA at these temperatures (see Figure [TJ. 
This time only the solution without a positivity constraint is 
shown since the regularization parameter for x 2 = 1 provides 
a positive solution. For comparison the forward fitted solutions 
from xrt_dem_iterative2 .pro (grey histograms) are shown 
and again match the regularized solution (and model DEM) over 
the majority of the temperature range but show a wider verti- 



cal spread at the smallest and largest temperatures. This consis- 
tent increase in error in the DEM solution from both methods 
suggests the poor response at these temperature extremes is the 
source of this uncertainty. The EM loci curves are also shown but 
again do not intersect at a single point as this is a broad DEM. 

In both examples shown so far only the zeroth-order con- 
straint has been used and in Figure [7] we show the resulting reg- 
ularized DEM for these cases for higher order constraint matri- 
ces. For the single Gaussian the higher order constraint removes 
the cluster of "noisy" data points at high temperature but also 
broadens the temperature resolution at around log 7q = 6.2-6.4. 
This shows that the constraint is changing the balance between 
these two temperature ranges and in this case the lower order 
constraint is preferable as it produces a better temperature reso- 
lution where the DEM is dominant. For the two Gaussian model 
DEM the higher order constraints again produce a broader tem- 
perature resolution which helps the regularized solution match 
the model better at low temperatures. Therefore in this case the 
higher order constraint is marginally preferable. As expected, 
in both cases the increase in order of the constraint matrix in- 
creases the "smoothing" of the regularized solution. In the sub- 
sequent analysis present here the zeroth order Lo constraint ma- 
trix is used throughout as it is generally sufficient to recover the 
expected DEM. 

We now consider the ability of the regularized inversion to 
recover single Gaussian model DEMs but with different widths 
and magnitudes, shown in Figure [8] Here we do not use the 
positivity constraint since for narrow DEMs the possible neg- 
ative regions are from temperatures where the contribution is 
tiny and consistent (within the errors) with zero and for broad 
DEMs the are no negative values recovered. As we set a = 1 the 
X 1 of the regularized solution is also approximately 1. Firstly, 
we increase the width of the Gaussian from the minimum ex- 
pected of crj - 0.1, the designed achievable temperature res- 
olution of SDO/AIA {Judge 2010), to cr T = 0.4, using a nor- 
malisation magnitude of Nq/( yfTno-j) = 10 23 cm _5 K _1 in Eq. 
fT2l shown in the top panels of Figure [8] The regularized in- 
version method is able to recover the model DEM at the limit 
of SDO/AIA's temperature resolution. For the widest model 
DEMs there is some slight deviation but it is consistent within 
the error bars. Again we compare the regularization method to 
xrt_dem_iterative2 .pro (blue histograms) and find similar 
solutions, with both methods deviating from the model DEM in 
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Fig. 8. DEMs recovered using regularized inversion (red error bars) and xrt_dem_iterative2 . pro (blue histograms, with 75 MC 
realisations) from simulated SDO/AIA data of Gaussian models (black dashed lines) of differing <x (increasing left to right) and 
normalisation magnitude (10 23 top row, 10 22 bottom row). 




6.0 6.5 7.0 6.0 6.5 7.0 6.0 6.5 7.0 6.0 6.5 7.0 

log,„T[MK] log 10 T[MK] log,„T[MK] log 10 T[MK] 



Fig. 9. DEMs recovered using regularized inversion (red error bars) and xrt_dem_iterative2 . pro (blue histograms, with 75 MC 
realisations) from simulated SDO/AIA data of models (black dashed lines) with a differing number of Gaussian components (left to 
right). 



similar temperature ranges, for instance below logT * 5.9 in the 
cr T = 0.4 model. We repeat the exercise but with model DEMs 
an order of magnitude smaller, No I '( y/27Tcr T ) — 10 22 cm" 5 K _1 , 
shown in the bottom panels of Figure [8] With the reduced signal 
to noise both methods produce poorer DEM solutions, partic- 
ularly for the narrowest Gaussian model showing the expected 
degradation of the resolution recoverable with increasing noise. 
In this case one of the xrt_dem_iterative2 .pro solutions de- 
viates greatly above logT w 6.7. For the broader DEMs the iter- 
ative forward fitting method produces DEMs closer to the model 
than the regularized inversion, though the later is mostly con- 
sistent within the indicated error bars. For the regularized inver- 
sion a better or "smoother" solution could be archived by using 
a higher order constraint matrix however with real data there is 
no prior knowledge to the form of the emission and therefore no 
indication if such a constraint is actually better. 

In Figure [9] we consider model DEMs constructed of sev- 
eral Gaussian components of width at the temperature reso- 
lution of SDO/AIA, <tj ~ 0.1. The first panel shows the 
same single Gaussian shown in Figure [8] for comparison. 
With two Gaussian components with centroid temperatures of 
logT = 6.0,7.0 (second panel in Figure [9]) the regularized 



solution closely matches the source model. The forward fit- 
ting approach (xrt_dem_iterative2 . pro) does get the general 
shape correct but produces a tall spread of solutions at the ex- 
tremes of the temperature range chosen and in between the two 
Gaussian components. The former will be due to the reduced 
SDO/AIA response in this range but for the discrepancy about 
logT = 6.5 this is not the case. The situation is even worse for 
xrt_dem_iterative2 .pro when three Gaussian components 
are considered (third panel in Figure [9] with centroid temper- 
atures of logT = 6.0, 6.5, 7.0) with the it producing a DEM so- 
lution completely different from the model. In one of the MC 
realisations the forward fitting method does recover two of the 
Gaussian components but still fails to recover the middle one. 
The regularized inversion method also has problems in recov- 
ering the model DEM but performs considerably better than 
xrt_dem_iterative2 .pro, producing a DEM with three dis- 
tinct components. Only between logT = 6.4-6.8 does it de- 
viate from the model, failing to match the second Gaussian 
peak and subsequent minima. In the final panel of Figure [9] 
we again consider three Gaussian components but this time 
with those of centroid temperatures log!T = 6.0,7.0 an order 
of magnitude smaller than the central component. This time 
xrt_dem_iterative2 . pro recovers the central component but 
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Fig. 10. DEMs recovered using regularized inversion (red error bars) from simulated Hinode/EIS data of Gaussian models (black 
dashed lines) of differing <x (increasing left to right). Here the error in the line intensity is taken as the Poisson noise, which range 
between 0.2 to 1 1% of the line intensity. In the bottom row the minimum of the EM loci curves (grey dotted lines) has been used as 
the guess solution %o(T). 
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Fig. 11. DEMs recovered using regularized inversion (red error bars) from simulated Hinode/EIS data of Gaussian models (black 
dashed lines) of differing <x (increasing left to right). Here the error is taken as 20% of the line intensity. In the bottom row the 
minimum of the EM loci curves (grey dotted lines) has been used as the guess solution %o(T). 



completely missing the two smaller ones. Again the regularized 
inversion performs far better recovering three distinct compo- 
nents though slightly overestimates the centroid position of the 
hottest component. 

3.2. Hinode/EIS Simulated Data 

Spectroscopic observations, such as those from the EUV 
Imaging Spectrometer (EIS) on Hinode, can potentially recover 
the DEM better than broadband multi-filter observations given 
the significant number of temperature sensitive lines available. 



The resulting observed line intensities and calculated contribu- 
tion functions can be easily used with our regularization method 
as it only requires an observable, associated error and response 
as a function of temperature for a number of filters or lines. Here 
we u se the atlas of EUV l ines (observable with Hinode/EIS) 
from lLandi & Youngl J2009) with 48 of them shown in the bot- 
tom panel of Figure Q] These are lines emitted at cooler tem- 
peratures than observed by SDO/AIA and so we will consider 
Gaussian models with centroid temperatures between logT = 
5.5 to 6.5. 
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Fig. 12. DEMs recovered using regularized inversion (red error bars) from simulated Hinode/EIS data of multiple component 
Gaussian models (black dashed lines) of differing cr (increasing left to right). Here the error in the line intensity is taken as the 
Poisson noise, which range between 0.2 to 1 1% of the line intensity. In the bottom row the minimum of the EM loci curves (grey 
dotted lines) has been used as the guess solution £o(T). 
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Fig. 13. DEMs recovered using regularized inversion (red error bars) from simulated Hinode/EIS data of multiple component 
Gaussian models (black dashed lines) of differing cr (increasing left to right). Here the error is taken as 20% of the line inten- 
sity. In the bottom row the minimum of the EM loci curves (grey dotted lines) has been used as the guess solution £o(T). 



The simulated line intensities are created as before with the 
pixel intensities for SDO/AIA but this time we consider not only 
the Poisson noise (6gi = -\fg) but also a systematic uncertainty 
that is a percentage of each line intensity. This corresponds to 
a temperature independent factor per line c,-, Sgi = gi/c/. We 
are using this approach as there are uncertainties in the relative 
(few percent), absolute (up to 22 %) and modelling (« 10%) of 
the contributio n functions for each l ine (private communication 
P. Young and lLandi & Y oung (2009)). Such a temperature inde- 
pendent systematic on the contribution function 5K, = Kj/c,, 
where c,- is the factor per line, is equivalent to the same system- 



atic on the observable 6gi = gi/ci, as discussed in ^2.51 So when 
calculating the error in the simulated observable for our chosen 
DEM we choose either the Poisson noise or a percentage error, 
in the latter case to represent this possibly dominant systematic 
uncertainty in the contribution functions. 

In Figure [10] we show the regularized inversion of Gaussian 
model DEMs of different widths using a Poisson uncertainty on 
the line intensity, again using the zeroth-order constraint and 
a — 1 (i.e. x 2 ~ !)■ F° r the model DEMs used the uncertain- 
ties given by 6gi = yfg range from 0.2 and 11% of the line 
intensities. The regularized solution very closely matches the 
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Fig. 14. Regularized DEMs (red error bars) and xrt_dem_iterative2 . pro solution (blue histograms, with 25 MC realisations) for 
simulated SDO/AIA data of CHIANTI DEM models (dashed black lines, left to right: quiet sun, active region, M-Class flare). 
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Fig. 15. Regularized DEMs (red error bars) for simulated Hinode/EIS lines of the CHIANTI quiet Sun DEM model. The error is 
taken as the Poisson noise and 1 % and 20% of the line intensity (left to right). In the bottom row the minimum of the EM loci curves 
(grey dotted lines) was used as the guess solution £o(T). 



model DEM for the narrowest cases (top left panels), match- 
ing the model down to the temperature resolution of the atomic 
data used to calculated the contribution functions, i.e. crj =0.1. 
For the wider Gaussian models (top right panels) the regularized 
solution recovers the model well about the peak emission but 
underestimates the emission at the lowest temperatures. These 
regularized DEMs were found using the two-stage approach of 
no initial guess solution %o(T) but we can also find the solutions 
using the minimum of the EM loci curves (grey dashed lines) as 
the guess solution (see 82. 2K shown in the bottom row of Figure 
[TUl For the narrowest Gaussian models there is little difference 
in the regularized solutions. The only improvement here is in the 
time it takes to perform the computation since the regularization 
has only been calculated once. However even with the two-stage 
approach the DEM is computed in just a few seconds. For the 
wider Gaussian models the regularized DEM found using the 
guess solution recovers the model DEM better, particularly at 
lower temperatures. 



This analysis is repeated but this time the error on the line 
intensities is taken to be 20% instead of the Poisson noise. With 
this increased uncertainty the regularization method does not re- 
cover the model DEMs as well, shown in Figure QT| For the 
narrowest DEM (top left panel) the regularized inversion recov- 
ers the majority of the model DEM but underestimates the peak 
emission and increases its width, the latter showing a reduction 
in temperature resolution due to noise. The use of the minimum 
of the EM loci curves as the guess solution %o(T) greatly im- 
proves the recovery of the source DEMs (bottom row) but still 
underestimates the peak emission in the narrowest cases. This 
demonstrates that the the temperature resolution had been inher- 
ently degraded by the increase in noise. 

We now consider model DEMs constructed of several 
Gaussian components of width at the temperature resolution of 
the atomic data, cr r w 0.1, shown with Poisson noise and then 
20% uncertainty in the line intensity in Figures[T2]and[T3] With 
the lower level of noise from the Poisson errors the regularized 
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Fig. 16. Regularized DEMs (red error bars) for simulated Hinode/EIS lines of the CHIANTI active region DEM model. The error 
is taken as the Poisson noise and 1% and 20% of the line intensity (left to right). In the bottom row the minimum of the EM loci 
curves (grey dotted lines) has been used as the guess solution %o(T). 



inversion recovers well the model DEMs (top row Figure ITSl. 
even in the cases with three Gaussian components. The use of 
the initial guess solution (bottom row Figure [T2l generally does 
worse in recovering the DEMs at the extremes of the tempera- 
ture ranges but produces similar results at the mid-range tem- 
peratures where the DEMs peak. With a larger error in the line 
intensity the regularized inversion struggles to recover the model 
DEMs (Figure \13i. In all cases it is able to recover the correct 
number of distinct components in the DEMs but is unable to 
match the peak emission, often producing flatter solutions. The 
use of the initial guess solution does help recover the model 
DEMs better (bottom row Figure IT3b but it still struggles with 
this noisy simulated data. 

4. Simulated Data: CHIANTI Model DEMs 

To test the regularization with more physically realistic DEMs 
we use those provided with the CHIANTI atomic database 
(iDere et al 1119971 120091) . We use t he DEMs for the qui et Sun 
(Idel Zannall 19991) . an a ctive region (lAndretta et al.ll2003l) and a 
M2-GOES Class flare (IDere &Cookll 19791) . With SDO/AIA we 
test all three of these models (see 34. Il l but for the Hinode/EIS 
lines we only consider the quiet Sun and active region models 
(see 34.2b since these lines are sensitive to temperatures predom- 
inantly below logT ss 6.5, lower than expected in a large flare. 
Again the regularized solutions were found using the zeroth- 
order constraint and as a — 1 the resulting solutions have ap- 
proximately^ 2 = 1. 



4.1. SDO/AIA Simulated Data 

The regularized DEMs recovered from simulated SDO/AIA data 
of the CHIANTI model DEMs is shown in Figure [H The first 
DEM shown is for the quiet Sun and the regularization recovers 
the DEM well but has very large error bars. This is understand- 
able given the faint, and hence noisy, emission: the peak emis- 



sion for the quiet Sun DEM is about 10 20 cm -5 K _1 which is 
over three orders of magnitude smaller that the Gaussian exam- 
ples in 33. II In comparison, xrt_dem_iterative2 .pro poorly 
recovers the model DEM producing a large spread of the MC re- 
alisations. The active region model DEM (middle panel Figure 
[T4T i is better recovered than the quiet Sun model which is ex- 
pected given that this DEM produces a stronger signal. Both 
methods have trouble in recovering the model DEM at the low- 
est temperatures (logT < 5.7) where the emission is small- 
est: the regularization solution underestimating the emission, 
the iterative approach overestimating. For the flare DEM (right 
panel Figure [T4T i both methods recover the model well about 
the peak temperature (logT * 7.0) will small vertical uncer- 
tainties. This is due to the large emission from the flare, about 
3 to 4 orders of magnitude larger than the quiet Sun and active 
region models. At lower temperatures the regularized solution 
matches the model DEM well though with a substantial hori- 
zontal spread. The xrt_dem_iterative2 .pro solutions devi- 
ate from the model below logr w 6.7 and produces a false peak 
about logT a; 6.1 when in the model it is a minimum. 

4.2. Hinode/EIS Simulated Data 

The regularized DEM of the CHIANTI quiet Sun model using 
the Hinode/EIS lines is shown in Figure [15] where it has been 
calculated using a variety of errors in the line intensities and 
with/without the minimum of the EM loci as the initial guess so- 
lution. When no initial guess solution is used (top row) the reg- 
ularization method is able to recover the main peak of emission 
(about \ogT « 6.0) for the three different error cases, Poisson 
noise, 1% and 20%. Only with the 1% errors on the line intensi- 
ties does the regularized solution properly recover that the emis- 
sion increases with decreasing temperature but does not match it 
very well (middle panel). Although it should be noted that this 
model produces very weak emission and the DEM recoverable 
from these multiple Hinode/EIS lines is considerably better than 
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Fig. 17. Regularized emissio n measure distributions (black solid line and grey error bars) for the active region core observations 
made by Warr en et al.l (HblO), their MCMC solution (red histogram) is shown for comparison. In the right-hand panels the positivity 
constraint was used and hence larger^ 2 values are achieved by the regularized solutions. The EM loci curves for the Hinode/EIS 
(coloured lines) and Hinode/XRT filter (black dashed line) are shown in the bottom left panel and their minimum was used as an 
initial guess solution for the regularizations shown in the bottom row. 



with the SDO/AIA simulated data (left panel Figure [14]). When 
the minimum of the EM loci curves is used as the initial guess 
solution (bottom row, Figure fT31 > the regularized solution does 
recover more of the model DEM but the problems with the low 
temperature component persist. 

In Figure[l6]we show the same analysis but this time for the 
CHIANTI active region DEM model. The 1% error case again 
well recovers the majority of model DEM but struggles at low 
and high temperatures when higher, and more realistic, errors 
in the line intensities are used. When the minimum of the EM 
loci curves are used as the initial guess solution (bottom row, 
FigurelToT) the regularization method does considerably better by 
impressively recovering the model DEM even in the situations 
with larger errors. 



5. Real Data: Hinode/XRT and EIS 

To test the performance of the regularization on real data we use 
the observ ations of an active r egion core with Hinode/EIS and 
XRT from lWarren et alJ (I2010I) . This article provides the line in- 
tensities, with errors (which were 20% to 2 3%), for 24 EIS line s 
and one XRT filter (Al-thick) (see Table 1, 1 Warren et ai](l2010l)) 
and co mputed the DEM (Figure 4, red histogram, IWarren et alJ 
(120101)') using the M CMC method packaged with PINTofALE 
(Kashyap & Drakdll998l) . Using the information given we use 
CHIANTI to calculate the contribution functions for the 24 EIS 
lines, and with the Hi node/XRT temper ature response function 
for the Al-thick filter dGolub et al.ll2007l) . and the quoted inten- 
sities and errors we calculate the zeroth-order regularized solu- 



tions and a = 1 (shown in Figure 1171). Here we have the regu- 
larized solutions found both with and without the positivity con- 
straint (left vs right columns) and the initial guess solution from 
the minimum of the EM loci curves (top vs bottom rows). Shown 
for co mparison is the MCMC solution found by IWarren et al.l 
d2010 *). Note that for all the regularized solutions the resulting 
R,jK matrix was almost diagonal over logT * 5.5-6.9 indicating 
that the regularization has successfully worked over this temper- 
ature range. With no positivity constraint (left-hand panels) the 
regularized solution is highly oscillatory between positive and 
negati ve values, an indication that this i s an over-regularized so- 
lution (ICraidll977t iBertero et alJ [T985). Although, the regular- 
ization does produce a maximum at the same temperature as the 
MCMC method. The use of the initial guess solution (bottom 
left panel Figure [T71 > results in minimal changes to the recov- 
ered DEM. The positivity constraint produces a closer match to 
the MCMC solution. When no initial guess solution is used (top 
right panel) a highly oscillatory DEM is recovered and deviates 
from the MCMC solution at the peak temperature and above. 
When the minimum of the EM loci curves are used as the ini- 
tial guess solution (bottom right panel) a smaller^ 2 is achieved 
while still having a positive solution. This regularized DEM is a 
close match to MCMC solution, even more so when one includes 
the sprea d of the DEM found from 250 MC solutions (shown in 
Figure 4, Warr en et al.l (1201 Oh ). The crucial difference though is 
in the computation time: the regularized solution is found within 
a few seconds, the MCMC method taking orders of magnitude 
longer. 
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6. Discussion & Conclusions 

In this work, we have applied a regularized inversion technique 
developed for RHESSI HXR analysis to multi-filter observations 
of hot solar plasma with Hinode/EIS, XRT and SDO/AIA. This 
method successfully recovers a variety of model DEMs from dif- 
ferent simulated broadband and spectroscopic data of varying 
noise and uncertainty. It is implementecQ using General Singular 
Value Decomposition (GSVD) and this has several advantages 
over previous a pproaches used to find th e regularized inversion 
of solar data (i.e lCraig|[T977l:lJudge et alJI 19971) : 

1 . It reliably recovers the DEM even from noisy data and us- 
ing only the weakest (yet mostly physically justifiable) con- 
straint, zeroth-order Lq; 

2. This method naturally determines the confidence interval for 
the regularized solution by calculating both the vertical and 
horizontal error bars allowing an objective assessment of 
the quality of the data, temperature response functions and 
DEM. The use of the R,iK matrix also provides additional 
information about the robustness of the regularized solution; 

3. It is computationally very quick, with DEMs recovered from 
SDO/AIA data in about a second or less and a few seconds 
when a large number of spectral lines are used, such as with 
Hinode/EIS. This is crucial considering the colossal amount 
of ever increasing data that is being accumulated by SDO, 
Hinode and expected from future missions. 

4. The recovery of a positive solution can be guaranteed though 
this approach can hide issues with the data and response 
functions. 

Misconceptions about the problems with inversion techniques 
(such as the smoothness criteria) still exist (e.g. lLandi et all 
1201 ll) but the examples shown in this paper clearly demonstrate 
that these views are misplaced and our regularized inversion ap- 
proach can robustly recover a variety of DEMs. Moreover the 
ability to easily determine the temperature uncertainty is a ma- 
jor improvement over other methods, crucial when trying to test 
the possibly isothermal nature of solar plasma. 

DEMs are a very useful tool for characterising the temper- 
ature distribution of the corona and trying to reveal the prop- 
erties of the mechanisms that are heating the plasma. However 
there are several caveats to the reconstructed DEMs which mean 
they should not be over-interpreted. Any method that recovers 
a DEM from a set of real data cannot guarantee that it is the 
actual solution given the possible uncertainties/errors in the data 
and temperature response functions (either instrumental or in the 
atomic physics). The use of more than one method is highly rec- 
ommended as this would highlight any issues and artifacts of a 
particular approach and demonstrate the role of the uncertainty 
in the data on the solution. The computational speed of our regu- 
larization method makes it a minor burden to include when using 
other techniques and also allows fast exploration of the effect of 
the uncertainties on the DEM, as demonstrated for the simulated 
data in this paper. Even if several techniques produce the same 
DEM, with similar uncertainties, this still leaves the possibil- 
ity of there being errors in the temperature response functions 
due to mis-calibration of the instrument or information in the 
atomic data being absent or erroneous. There is a continuous 
effort to improve the atomic physics data, through the sterling 
work of the CHIANTI team and others. At the moment there 
are known issues with the SDO/AIA response functions (e.g. 

3 The code written in IDL requiring SSW is available online: 
http://www.astro.gla.ac.uk/~iain/demreg/ 



lAschwanden & Boernerl201 lh and the effect of this can be tested 
by simulating SDO/AIA data using the "correct" response but 
recover the DEM from this data using the erroneous responses. 
However if the emission is not optically thin or not in thermal 
equilibrium the DEM and response functions are not appropri- 
ate to describe the temperature distribution of the plasma. 

All these problems are compounded by the vast quantities 
of solar data now available resulting in the current need for a 
variety of tools to investigate DEMs quickly, transparently and 
easily. The use of simulated data is a vital approach to investi- 
gate the role of each of these issues and combined with a com- 
putationally fast algorithm that provides error estimates, such as 
our regularization method, should provide a quicker determina- 
tion of the reliability of DEMs and what can be interpreted from 
them. 

Acknowledgements. This work is supported by a STFC grant (IGH, EPK). 
Financial support by the European Commission through the FP7 HESPE 
Network is gratefully acknowledged. We would like to thank the referee for their 
constructive criticisms that helped us greatly improve this paper. 



References 

Andretta, V., Del Zanna, G., & Jordan, S. D. 2003, A&A, 400. 737 
Aschwanden, M. J. & Boemer, P. 201 1, ApJ, 732, 81 
Bertero, M., DeMol, C, & Pike, E. R. 1985, Inverse Problems, 1, 301 
Bertero, M., DeMol, C, & Pike, E. R. 1988, Inverse Problems, 4, 573 
Bertero, M. & Dovi, V. 1981, Optica Acta, 28, 1635 
Boemer, P., Edwards, C, Lemen, J., et al. 2011, Sol. Phys., 193 
Brosius, J. W., Davila, J. M., Thomas, R. J., & Monsignori-Fossi, B. C. 1996, 
ApJS, 106, 143 

Brown, J. C, Emslie, A. G., Holman, G. D., et al. 2006, ApJ, 643, 523 
Craig, I. J. D. 1977. A&A, 61, 575 

Craig, I. J. D. & Brown, J. C. 1986, Inverse problems in astronomy (Bristol: 
Hilger) 

Culhane, J. L., Harra, L. K., James, A. M., et al. 2007, Sol. Phys., 243, 19 
de Villiers, G. D., McNally, B., & Pike, E. R. 1999, Inverse Problems, 15, 615 
del Zanna, G. 1999, PhD thesis, , Univ. of Central Lancashire, (1999) 
Dere, K. P. & Cook, J. W. 1979, ApJ, 229, 772 

Dere, K. P., Landi, E., Mason, H. E., Monsignori Fossi, B. C, & Young, P. R. 

1997, A&AS, 125, 149 
Dere, K. P., Landi, E., Young, P. R., et al. 2009, A&A, 498, 915 
Fletcher, L., Dennis, B. R., Hudson, H. S., et al. 2011, Space Sci. Rev., 159, 19 
Golub, L., Deluca, E., Austin, G., et al. 2007, Sol. Phys., 243, 63 
Golub, L., Deluca, E. E., Sette, A., & Weber, M. 2004, in Astronomical Society 

of the Pacific Conference Series, Vol. 325, The Solar-B Mission and the 

Forefront of Solar Physics, ed. T. Sakurai & T. Sekii, 217 — l- 
Goryaev, F. F, Parenti, S., Urnov, A. M., et al. 2010, A&A, 523, A44+ 
Hannah, I. G., Hudson, H. S.. Battaglia, M., et al. 2011, Space Sci. Rev., 159, 

263 

Hansen, P. C. 1992, Inverse Problems, 8, 849 
Hirayama, T. 1974, Sol. Phys., 34, 323 
Judge, P. G. 2010, ApJ, 708, 1238 

Judge, P. G., Hubeny, V., & Brown, J. C. 1997, ApJ, 475, 275 
Kashyap, V. & Drake, J. J. 1998, ApJ, 503, 450 

Kontar, E. P., Piana, M., Massone, A. M., Emslie, A. G., & Brown, J. C. 2004, 

Sol. Phys., 225, 293 
Landi, E., Reale, F, & Testa, P. 201 1, ArXiv e-prints 
Landi, E. & Young, P. R. 2009, ApJ, 706, 1 

Lemen, J. R., Title, A. M., Akin, D. J., et al. 201 1, Sol. Phys., 172 
Lin. R. P., Dennis, B. R., Hurford, G. J., et al. 2002, Sol. Phys., 210, 3 
Mariska, J. T. 1992, The solar transition region, ed. Mariska, J. T. 
Massone, A. M., Emslie, A. G., Kontar, E. P., et al. 2004, ApJ, 613, 1233 
Mcintosh, S. W., Charbonneau, P., & Brown, J. C. 2000, ApJ, 529, 1115 
Monsignori Fossi, B. C. & Landini, M. 1991, Advances in Space Research, 11, 
281 

Monsignori-Fossi, B. C. & Landini, M. 1992, Mem. Soc. Astron. Italiana, 63, 
767 

Morozov, V. 1967, Doklady Akademii Nauk SSSR, 175, 1225 
O'Dwyer, B., Del Zanna. G., Mason, H. E., et al. 201 1, A&A, 525, A137+ 
Parenti, S., Bromage, B. J. I., Poletto, G., et al. 2000, A&A, 363, 800 
Parker, E. N. 1988, ApJ, 330, 474 



14 



Hannah & Kontar: Regularized Inversion of Hinode and SDO data 



Phillips, K. J. H., Feldman, U., & Landi, E. 2008, Ultraviolet and X-ray 

Spectroscopy of the Solar Atmosphere, ed. Phillips, K. J. H., Feldman, U., 

& Landi, E. (Cambridge University Press) 
Piana, M. & Bertero, M. 1997, Inverse Problems, 13, 441 
Piana, M., Massone, A. M., Kontar, E. P., et al. 2003, ApJ, 595, L127 
Prato, M., Piana, M., Brown, J. C, et al. 2006, Sol. Phys., 237, 61 
Press, W. H., Teukolsky, S. A., Vetterling, W. T., & Flannery, B. P. 1992, 

Numerical recipes in C. The art of scientific computing, ed. Press, W. H., 

Teukolsky, S. A., Vetterling, W. T., & Flannery, B. P. 
Reale, F., Testa, P., Klimchuk, J. A., & Parenti, S. 2009, ApJ, 698, 756 
Schmelz, J. T., Jenkins, B. S., Worley, B. T., et al. 201 1, ApJ, 731, 49 
Schmitt, J. H. M. M., Drake, J. J., Stern, R. A., & Haisch, B. M. 1996, ApJ, 457, 

882 

Siarkowski, M., Falewicz, R., Kepa, A., & Rudawy, P. 2008, Annales 

Geophysicae, 26, 2999 
Tikhonov, A. N. 1963, Soviet Math. Dokl., 4, 1035 
Warren, H. P., Winebarger, A. R., & Brooks, D. H. 2010, ApJ, 71 1, 228 
Weber, M. 2009, in Astronomical Society of the Pacific Conference Series, 

Vol. 415, The Second Hinode Science Meeting: Beyond Discovery-Toward 

Understanding, ed. B. Lites, M. Cheung, T. Magara, J. Mariska, & K. Reeves, 

32-+ 

Weber, M. A., Deluca, E. E., Golub, L., & Sette, A. L. 2004, in IAU 
Symposium, Vol. 223, Multi- Wavelength Investigations of Solar Activity, ed. 
A. V. Stepanov, E. E. Benevolenskaya, & A. G. Kosovichev, 321-328 

Weber, M. A., Schmelz, J. T., DeLuca, E. E., & Roames, J. K. 2005, ApJ, 635, 
L101 

Woods, T. N., Eparvier, F. G., Hock, R., et al. 2010, Sol. Phys., 3 



